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Abstract 



J2' 



We consider direct diphoton production in hadron collisions. We compute the next-to-next-to-leading order (NNLO) 
QCD radiative corrections at the fully-differential level. Our calculation uses the qj subtraction formalism and it is 
implemented in a parton level Monte Carlo program, which allows the user to apply arbitrary kinematical cuts on the 
final-state photons and the associated jet activity, and to compute the corresponding distributions in the form of bin 
histograms. We present selected numerical results related to Higgs boson searches at the LHC and the Tevatron, and 
we show how the NNLO corrections to diphoton production are relevant to understand the main background of the 
decay channel H — > yy of the Higgs boson H. 
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1. Introduction 

Diphoton production is a relevant process in hadron col- 
lider physics. It is both a classical signal within the 
Standard Model (SM) and an important background for 
Higgs boson and new-physics searches. The origin of 
the Electroweak symmetry breaking is currently being 
investigated at the LHC by searching for the Higgs bo- 
son and eventually studying its properties. Recent re- 
sults in the search for the SM Higgs Boson at the LHC 
indicates the observation of a new particle iQISL which 
is a neutral boson with a mass M ~ 125 GeV. In 
this spectacular new observation, as well as in previous 
searches and studies, the preferred search mode involves 
Higgs boson production via gluon fusion followed by 
the rare decay into a pair of photons. Therefore, it is 
essential to count on an accurate theoretical description 
of the various kinematical distributions associated to the 
production of pairs of prompt photons with large invari- 
ant mass. 

We are interested in the process pp — > yyX, which, at 
the lowest order in perturbative QCD, occurs via the 
quark annihilation subprocess qq — > yy. The QCD 
corrections at the next-to-leading order (NLO) in the 
strong coupling as involve the quark annihilation chan- 
nel and a new partonic channel, via the subprocess 
18 ~ > YY a - These corrections have been computed and 
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implemented in the fully-differential Monte Carlo codes 
DIPHOX H, 2gammaMC [% and MCFM |H]. A calcula- 
tion that includes the effects of transverse-momentum 
resummation is implemented in RESBOS |0]. At the 
next-to-next-to-leading order (NNLO), the gg channel 
starts to contribute, and the large gluon-gluon luminos- 
ity makes this channel sizeable. Part of the contribution 
from this channel, the so called box contribution, was 
computed long ago [7] and its size turns out to be com- 
parable to the lowest-order result. 

Besides their direct production from the hard subpro- 
cess, photons can also arise from fragmentation sub- 
processes of QCD partons. The computation of frag- 
mentation subprocesses requires (poorly known) non- 
perturbative information, in the form of parton frag- 
mentation functions of the photon. The complete NLO 
single- and double-fragmentation contributions are im- 
plemented in DIPHDX [3]. The effect of the fragmenta- 
tion contributions is sizeably reduced by the photon iso- 
lation criteria that are necessarily applied in hadron col- 
lider experiments to suppress the very large irreducible 
background (e.g., photons that are faked by jets or pro- 
duced by hadron decays). The standard cone isolation 
and the 'smooth' cone isolation proposed by Frixione 
0] are two of these criteria. The standard cone isola- 
tion is easily implemented in experiments, but it only 
suppresses a fraction of the fragmentation contribution. 
The smooth cone isolation (formally) eliminates the en- 
tire fragmentation contribution, but its experimental im- 
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plementation is still under study |@] . However, it is im- 
portant to anticipate (work to appear) that in some kine- 
matical regions (e.g for Higgs boson searches) QCD cal- 
culations that apply the standard cone and the Frixione 
isolation criteria give basically very similar quantitative 
result^. 



2. Diphoton production at NNLO 

We consider the inclusive hard-scattering reaction 



hi + hi 



77 



X 



(1) 



where the collision of the two hadrons, hi and hz, pro- 
duces the diphoton system F = yy with high invariant 
mass M yj . The evaluation of the NNLO corrections to 
this process requires the knowledge of the correspond- 
ing partonic scattering amplitudes with X = 2 partons 
(at the tree level lllOll ). X = 1 parton (up to the one- 
loop level [11]) and no additional parton (up to the two- 
loop level 11211 ) in the final state. The implementation of 
the separate scattering amplitudes in a complete NNLO 
(numerical) calculation is severely complicated by the 
presence of infrared (IR) divergences that occur at in- 
termediate stages. The qj subtraction formalism 11311 
is a method that handles and cancels these unphysical 
IR divergences up to the NNLO. The formalism applies 
to generic hadron collision processes that involve hard- 
scattering production of a colourless high-mass system 
F. Within that framework [13], the corresponding cross 
section is written as: 



[da 

d(T' 



F+jets 
(N)LO 
CT 

(N)LO 



] , (2) 



where dcr^^ represents the cross section for the pro- 
duction of the system F plus jets at (N)LO accuracy^, 
and dcr^ )L0 is a (IR subtraction) counterterm whose ex- 
plicit expression lfl5ll is obtained from the resumma- 
tion program of the logarithmically-enhanced contri- 
butions to qj distributions. The 'coefficient' H F N ^ NL0 , 
which also compensates for the subtraction of dcr^ LO , 
corresponds to the (N)NLO truncation of the process- 
dependent perturbative function 
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(3) 



1 The use of the same parameters in both criteria is understood. 
2 In the case of diphoton production, the NLO calculation of 
dcrJJ^Q^ was performed in Ref. I4|. 



The NLO calculation of dcr F requires the knowledge of 
f H Fm , and the NNLO calculation also requires f H F{1 \ 
The general structure of <H FW is explicitly known 01611 : 
exploiting the explicit results of <H Ft - x> for Higgs Il3lll7ll 
and vector boson fl8ll production, we have generalized 
the process-independent relation of Ref. [16] to the cal- 
culation of the NNLO coefficient <H F(Y > '. 



3. Quantitative results 

We have performed our fully-differential NNLO calcu- 
lation 1 19] of diphoton production according to Eq. ©. 
The NNLO computation is encoded in a parton level 
Monte Carlo program, in which we can implement ar- 
bitrary IR safe cuts on the final-state photons and the 
associated jet activity. We concentrate on the direct pro- 
duction of diphotons, and we rely on the smooth cone 
isolation criterion (8J[. Considering a cone of radius 
r - ^{Arf) 2 + (A<p) 2 around each photon, we require 
the total amount of hadronic (partonic) transverse en- 
ergy Ej inside the cone to be smaller than E^L(f), 



E T <E<f max (r) = E Tmax 



1 - cos r 
1 - cos R 



(4) 



where Ej max can be a fixed value or a fraction of the 
transverse momentum of the photon (Ej max - £y p y T , 
with < € y < 1); the isolation criterion Ej < E'j^Jj) 
has to be fulfilled for all cones with r < R. We use the 
MSTW 2008 S sets of parton distributions, with den- 
sities and as evaluated at each corresponding order, and 
we consider Nf = 5 massless quarks/antiquarks and glu- 
ons in the initial state. The default renormalization (jUr) 
and factorization [fif ) scales are set to the value of the 
invariant mass of the diphoton system, /jr = [if — M yy . 
The QED coupling constant a is fixed to a — 1/137. 

To present some quantitative results, we consider dipho- 
ton production at the LHC ( -y^ = 14 TeV). We apply 
typical kinematical cuts used by ATLAS and CMS Col- 
laborations in their Higgs boson search studies. We re- 
quire the harder and the softer photon to have transverse 
momenta p^ 1 > 40 GeV and p s ° fla > 25 GeV, re- 
spectively. The rapidity of both photons is restricted to 
|y y | < 2.5, and the invariant mass of the diphoton sys- 
tem is constrained to lie in the range 20 GeV < M yy < 
250 GeV. The isolation parameters are set to the val- 
ues e y = 0.5, n — 1 and R = 0.4. Performing the 
QCD computation, we observe lfl9ll that the value of 
the cross section remarkably increases with the pertur- 
bative order of the calculation. This increase is mostly 
due to the use of very asymmetric (unbalanced) cuts on 
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Figure 1 : Invariant mass distribution of the photon pair at the LHC = 
14 TeV): LO (dots), NLO (dashes) and NNLO (solid) results. We also present 
the results of the box and NLO+box contributions. The inset plot shows the 
corresponding K-factors. 



Figure 2: Diphoton cross section as a function of the azimuthal separation 
of the two photons. Data from CMS |21] ( = 7 TeV) are compared to the 
NNLO calculation Q3]. 



the photon transverse momenta. At the LO, kinematics 
implies that the two photons are produced with equal 
transverse momentum and, thus, both photons should 
have p y T > 40 GeV. At higher orders, the final-state ra- 
diation of additional partons opens a new region of the 
phase space, where 40 GeV > p s ° ftel > 25 GeV. Since 
photons can copiously be produced with small trans- 
verse momentum 01911 . the cross section receives a size- 
able contribution from the enlarged phase space region. 
This effect is further enhanced by the opening of a new 
large-luminosity partonic channel at each subsequent 
perturbative order. In Fig. Q] we show the LO, NLO 
and NNLO invariant mass distributions at the default 
scales. We also plot the gluonic box contribution (com- 
puted with NNLO parton distributions) and its sum with 
the full NLO result. The inset plot shows the K-factors 
defined as the ratio of the cross sections at two subse- 
quent perturbative orders. We note that }<^nw/NW j s 
sensibly smaller than K NLO/LO , and this fact indicates 
an improvement in the convergence of the perturbative 
expansion. We find that about 30% of the NNLO cor- 
rections is due to the gg channel (the box contribution 
is responsible for more than half of it), while almost 
60% still arises from the next-order corrections to the 
qg channel. The NNLO calculation includes the pertur- 
bative corrections from the entire phase space region (in 
particular, the next-order correction to the dominant qg 
channel) and the contributions from all possible partonic 
channels (in particular, a fully-consistent treatment of 
the box contribution to the gg channel 0). Owing to 



these reasons, the NNLO result can be considered a re- 
liable estimate of direct diphoton production, although 
further studies (including independent variations of /i« 
and fif , and detailed analyses of kinematical distribu- 
tions) are necessary to quantify the NNLO theoretical 
uncertainty. 



Recent results from the LHC 112 ll 12211 and the Teva- 
tron i23ll show some discrepancies between the data 
and NLO theoretical calculations of diphoton produc- 
tion. Basically, discrepancies were found in kinematical 
regions where the NLO calculation is effectively a LO 
theoretical description of the process. Such phase space 
regions (away from the back-to-back configuration) are 
accessible at NLO for the first time, due to the final-state 
radiation of the additional parton 0. Figure [2] shows a 
measurement by CMS 02 ill , of the diphoton cross sec- 
tion as a function of the azimuthal angle A(p yy between 
the photons. The data are compared with our NLO and 
NNLO calculations lfl9ll . The acceptance criteria used 
in this analysis ( ys = 7 TeV) require: p^ del > 23 GeV 
and /7^? fter > 20 GeV. The rapidity of both photons is 
restricted to \y y \ < 2.5, and the invariant mass of the 
diphoton system is constrained to be M yy > 80 GeV. 
The isolation parameters have the values e y = 0.05, 
n = 1 and R = 0.4. 

The histograms in Fig. [2] show that the NNLO QCD re- 
sults remarkably improve the theoretical description of 
the CMS data throughout the entire range of A(f> yy . 

In Figure [3] we present the invariant mass distribu- 
tion for diphoton production at the Tevatron (^Js — 



3 The calculation | 4] of the next-order gluonic corrections to the 
box contribution indicates an increase of the NNLO result by less than 
10% if Myy > 100 GeV. 



The low-mass region (M yy 
this case. 



SOGeV) in Figure[T]also belongs to 
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1 .96 TeV) calculated at NNLO, compared with a mea- 
surement performed by CDF l23ll . The acceptance 
criteria in this case require: pj dld " > 17 GeV and 
^softer > 25 GeV. The rapidity of both photons is re- 
stricted to \y 7 \ < 1. The isolation parameters are set 
to the values Ej max = 2 GeV, n — 1 and R = 0.4, 
and the minimum angular separation between the two 



photons is R 



yy 



0.4. Though in this case the in- 
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Figure 3: Diphoton cross section as a function of the invariant mass of the two 
photons. Data from CDF E3l (V? = 1.96 TeV) are compared to the NNLO 
calculation. 



crease from the LO to the NLO result is considerably 
smaller than at the LHC HI, the NNLO QCD correc- 
tions still improve remarkably the theoretical descrip- 
tion of the CDF data, in particular in the low mass re- 
gion (M yy < 2p'™ rder = 34 GeV). 

We note that the CMS and CDF data are obtained by 
using the standard cone isolation criterion and the con- 
straint in Eq. (0 is applied only to the cone of radius 
r — R. Since the smooth isolation criterion used in our 
calculation (we apply Eq. (0| for all cones with r < R) 
is stronger than the photon isolation used by CMS and 
CDF, we remark that our NLO and NNLO results can- 
not overestimate the corresponding theoretical results 
for the experimental isolation criterion. 
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The results illustrated in this contribution show that the 
NNLO description of diphoton production is essential 
to understand the phenomenology associated to this pro- 
cess, and therefore, the NNLO calculation is a relevant 
tool to describe the main background for Higgs boson 
searches and studies. 
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